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We present an approach to characterize genuine multiparticle entanglement using appropriate 
approximations in the space of quantum states. This leads to a criterion for entanglement which 
can easily be calculated using semidefinite programming and improves all existing approaches sig- 
nificantly. Experimentally, it can also be evaluated when only some observables are measured. 
Furthermore, it results in a computable entanglement monotone for genuine multiparticle entangle- 
ment. Based on this, we develop an analytical approach for the entanglement detection in cluster 
states, leading to an exponential improvement compared with existing schemes. 
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Introduction — The characterization of multiparti- 
cle quantum correlations is relevant for many physical 
systems like atoms in optical lattices, superconducting 
qubits or nitrogen-vacancy centers in diamond, to only 
name some recent examples In the field of quantum 
information, multiparticle entanglement is viewed as a re- 
source, enabling tasks like measurement-based quantum 
computation U or high-precision metrology Q • In spite 
of many efforts, the characterization of these correlations 
turns out to be difficult. Especially genuine multipartite 
entanglement, which is most important from the experi- 
mental point of view, remains unruly and only scattered 
results concerning its characterization are known • 

In this Letter, we derive a general method to char- 
acterize genuine multiparticle entanglement using suit- 
able relaxations. This relaxed problem turns out to be 
good-natured, can be tackled with different methods and 
results in a criterion that can be considered as a general- 
ization of the Peres-Horodecki criterion || to the multi- 
partite case. The goal of our work is two-fold: First, 
we present powerful criteria for genuine multiparticle 
entanglement, which can be efficiently evaluated using 
semidefinite programming and improve existing condi- 
tions significantly. They work for multi-qubit systems, 
continuous-variable or hybrid systems and can be evalu- 
ated, even if the mean values of only a few observables are 
known. Furthcmore, they lead to a computable entangle- 
ment monotone for genuine multiparticle entanglement. 

Second, our method allows to analytically derive en- 
tanglement conditions for the family of cluster states , 
which arc important states for tasks like measurement- 
based quantum computation. The sensitivity of these 
conditions improves exponentially with the number of 
qubits, which is an exponential gain compared with the 
existing conditions. As a sideproduct of our investiga- 
tions, we will also estimate the volume of the set of gen- 
uinely multipartite entangled states and gain insight into 
the geometrical form of the set of biseparable states. 

Situation — We start by considering a three-particle 




FIG. 1: For three qubits, there are three convex sets of states 
that are separable with respect to a fixed bipartition, namely 
the bipartitions A\BC, B\AC and C\AB (blue, dashed lines). 
Their convex hull (thick blue, dashed line) is the set of bisep- 
arable states. Each of the three sets is contained within the 
set of states that are PPT with respect to the corresponding 
bipartition (red, solid lines). Their convex hull forms the set 
of PPT mixtures (thick red, solid line). 

quantum state g. This state is separable with re- 
spect to some bipartition, say, A\BC, if it is a mix- 
ture of product states with respect to this bipartition, 

8 = E fc 9fcl^)('/ , Al ® H>bg){^bc\^ where thc 9k form a 
probability distribution. We denote these states by 
@A\bc Similarly, we can define the separable states for 
the two other possible bipartitions, Q S ^ AC and Q^ AB - 

Then, a state is called biseparable if it can be written 
as a mixture of states which are separable with respect 
to different bipartitions Q. That is, one has 

bs scp , scp . sep / -i \ 

Q =P10 A IBC+P^B\AC + P^C\AB- W 

On the other hand, a state that is not biseparable is 
called genuinely multipartite entangled. Whenever we 
talk about multipartite entangled states in the following, 
we refer to genuinely multipartite entangled states. 
To characterize multipartite entanglement, we apply 
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the method illustrated by Fig. [TJ Instead of states like 
@A\bc that are separable with respect to a fixed bipar- 
tition, we consider a larger set of states, which can be 
more easily characterized. For instance, for the biparti- 
tion A\BC we may consider states which have a positive 
partial transpose (PPT) [lOfl . It is well-known that sep- 
arable states are also PPT [8j . We denote such states by 
S V A\bc ( an d analogously for the other bipartitions) . 
Thus, we ask whether a state can be written as 



pmix _ ppt 

Q — PiQa\bc 



P2Q P B\AC 



P30 V C\AB- ( 2 ) 



We call states of this form PPT mixtures. Clearly, any 
biscparable state is a PPT mixture, so proving that a 
state is no PPT mixture implies genuine multipartite en- 
tanglement. There are examples of states, which arc PPT 
with respect to any bipartition, but nevertheless multi- 
partite entangled [111 ]. Hence, not all multipartite entan- 
gled states can be detected in this way, but, as we will 
see, often the set of PPT mixtures is a very good approxi- 
mation to the set of biseparable states. Finally, note that 
all definitions can be extended to TV particles. Also, one 
may use other relaxations of bipartite separability, e.g. 
apply the criterion of Doherty et al. [l2| . 

The advantage of considering PPT mixtures instead of 
biseparable states is that the set of PPT mixtures can 
be fully characterized with the method of linear semidef- 
inite programming (SDP) (l3j — a standard problem of 
constrained convex optimization theory. Moreover, PPT 
mixtures can also be characterized analytically. 

Characterization via entanglement witnesses — An 
entanglement witness is an observable W that is non- 
negative on all biseparable states, but has a negative 
expectation value on at least one entangled state. Let 
us first consider two particles, A and B. Then a decom- 
posable witness is a witness W that can be written as 
W = P + Q Ta , where P and Q have no negative eigen- 
values, (they are positive semidefinite, P, Q > 0) and Ta 
is the partial transpose with respect to A |l4| . 

For more than two particles, we call a witness W fully 
decomposable if, for every subset M of all systems, it is 
decomposable with respect to the bipartition given by M 
and its complement M. This means, there exist positive 
semidefinite operators Pm and Qm, such that 



for all M : W = P L 



M 



Qm 



(3) 



This observable is positive on all PPT mixtures, as it is 
non-negative on all states which are PPT with respect to 
some bipartition. But also the converse holds: 

Observation. If g is not a PPT mixture, then there 
exists a fully decomposable witness W that detects g. 

Proof — The set of PPT mixtures is convex and com- 
pact. Therefore, for any state outside of it, there is a wit- 
ness that is positive on all PPT mixtures. Furthermore, 
positivity on all states that arc PPT with respect to a 
fixed (but arbitrary) bipartition implies that the witness 



is decomposable with respect to this fixed (but arbitrary) 
bipartition [H]. Thus, W = Pm + Q T A f for any M. □ 
Practical evaluation — To find a fully decomposable 
witness for a given state, the convex optimization tech- 
nique SDP becomes important, since it allows us to op- 
timize over all fully decomposable witnesses. Given a 
multipartite state g, the search is given by 

min Tt(Wq) (4) 
s.t. Tr(W) = 1 and for all M : 

W = P m + Qm*,Qm>0, Pm>0. 

The free parameters are given by W and the operators 
Pm for every subset M. If the minimum in Eq. @ is 
negative, g is not a PPT mixture and hence genuinely 
multipartite entangled. The operator W for which the 
negative minimum is obtained is a fully decomposable 
witness. Note that, due to X ThI = (X T ) T ~ and X > 
X T > 0, a witness that is decomposable with respect to 
M is also decomposable with respect to M. Thus, one 
needs to check only half of all subsets in practice. 

Eq. Q has the form of a semidefinite program fl3j |. 
In contrast to usual optimization problems, global opti- 
mality of an SDP can be certified and the solution can 
efficiently be computed via interior-point methods. In 
practice, Eq. Q can be solved with few lines of code, us- 
ing e.g. the parser YALMIP [l6| and, as solvers, SeDuMi 



17] or SDPT3 18]. Our implementation in MATLAB 



named PPTMixer can be found online [15| . 

Let us discuss two variations of Eq. (U). First, in order 
to reduce the number of free parameters, one can restrict 
oneself to witnesses W that obey W Tm > for all M, 
i.e., Pm = for all M. In the following, we will call 
these witnesses fully PPT witnesses. For bipartite sys- 
tems, decomposable witnesses and fully PPT witnesses 
detect the same states. For the multipartite case, fully 
PPT witnesses arc not as good as fully decomposable 
witnesses, but they are easier to characterize. 

Second, this SDP can also be modified to account 
for the case that, instead of a full tomography, only 
a restricted set of observables is measured. Let O = 
{Oi,...Ok} be such a set of observables. Then, adding 
W = Yji=i to the constraints in Eq. (01 results in 
an SDP that searches for witnesses which can be evalu- 
ated knowing these observables. Note that for this pro- 
gram the free parameters are given by the real numbers Xt 
and their number might be considerably smaller than in 
Eq. (|4]). If the minimum in Eq. (j4|) is non- negative, there 
exists a PPT mixture with expectation values (Oj). How- 
ever, one may then add further observables to O and run 
the SDP again. Repeating this procedure gives more and 
more sensitive tests. We will discuss an example later. 
In practice, this program can even decide separability if 
the Oi have already been measured, so it can be used to 
gain new insights into already performed experiments. 
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state 


white noise tolerances ptoi 




fully decomposable 


before 


\GHZ 3 )* 


0.571 


0.571 [JJ 


\GHZ A )* 


0.533 


0.533 [JJ 


\w 3 y 


0.521 


0.421 [JJ 


\Wa) 


0.526 


0.444 [JJ 


\cuy 


0.615 


0.533 [25J 


\D2A) 


0.539 


0.471 [26] 


l*S,4> 


0.553 


0.317 [2JJ 



TABLE I: White noise tolerances of the fully decomposable 
witnesses obtained by the SDP of Eq. Q compared with the 
corresponding tolerances of the most robust criteria known 
so far. For the states marked by *, we verified that adding 
more white noise than what is tolerated by Eq. Q results in 
a biseparable state, so the values are optimal. 



But before proceeding to the examples, let us note 
three more facts. First, in the formulation no dimension 
of the Hilbert space is fixed. Consequently, our approach 
is valid for any dimension and combined with the meth- 
ods of Ref. [19J it can be directly used to study multipar- 
tite entanglement in continuous-variable or hybrid sys- 
tems poj . For continuous variables, it can be employed 
complementary to the methods of Ref. [2l[ . 

Second, our approach can also be used to quantify gen- 
uine multipartite entanglement. If in Eq. Q the trace 
normalization Tr(W) = 1 is replaced by < Pm < 1 
and < Qm < 1, the negative witness expectation value 
is a multipartite entanglement monotone, since it obeys 
the following properties: (i) It vanishes on all biseparable 
states, (ii) It is convex, (iii) The quantity does not in- 
crease under protocols that consist of local operations of 
each party and classical communication between them, 
(iv) It is invariant under local basis changes. While most 
of these properties are straightforward to see — in par- 
ticular, (iv) is implied by (iii) — , the proof of property 
(iii) is more technical (cf. Appendix A). Note that, in 
the bipartite case, this monotone becomes the negativ- 
ity 

Third, as mentioned before, there are other possible 
choices of supersets for the set of separable states, e.g. 
the set of states that have a symmetric extension on a 
larger Hilbert space (cf. Appendix B and fl2j]). 

Numerical examples — We test the criterion of Eq. (|4]) 
for important pure three- and four-qubit states prepared 
in many experiments 23[, using the white noise tol- 
erance as a figure of merit. It is given by the max- 



imal amount p to \ of white noise for which the state 
fbtoi) = (1 - Pto\)\i>){ip\ +ptoil/2 n is still detected as 
entangled [24|. Table Q] shows the white noise tolerances 
of our criterion, compared with the most robust criteria 
so far. 

Strikingly, the tolerances of the witnesses obtained by 
our SDP arc significantly higher than previous ones. For 



the GHZ and the W state of three qubits and the GHZ 
and the linear cluster state of four qubits, we even obtain 
the best white noise tolerance possible, since we are able 
to show that for a larger amount of white noise the state 
becomes biseparable (see Appendix C). This shows that 
our criterion is indeed optimal for these cases. 

To show that the criterion of Eq. Q works well 
for a restricted set of obscrvablcs, we consider the 
four-qubit Dicke state with two excitations |r>2,4) [Hj]. 
For this state, the SDP yields a witness Wd ( scc 
Appendix D) that consists of the observables O = 
{X®\ Y®\ Z®\ X X X 2 Y^ X X X % Z^ Y 1 Y 2 Z 3 Z 4 }, 
their distinct permutations and other observables that 
can be measured in the same run. For example, a 
local measurement of ATiA^A^A^ yields knowledge of 
the expectation value of X1X2IL3X4. The SDP finds 
a witness consisting of 0\ = X® 4 , O2 = Y® A and 
observables obtained by replacing some Paulis by the 

identity. Already with these observables, the white noise 

(2) 

tolerance is p y to { rj 0.29495. We can proceed in this way 
and use additional observables Oi from the set O - 
including their permutations and observables obtained 
by replacing Pauli operators by 1 — to produce strictly 
1 Their white noise tolerances p[ 1 ^ 
0.45200 and 



stronger witnesses W D . 
arc p[ 3 J « 0.38379, p£> 



finally pf^ 



0.38383, 



0.53914 as in Tabic H since W D = W, 



(6) 



Third, we compute a lower bound on the volume of 
genuinely multipartite entangled states. We created sam- 
ples of 10 4 random mixed three-qubit states uniformly 
distributed with respect to the Hilbert-Schmidt distance 
(or the Bures distance) and check whether they are gen- 
uinely multipartite entangled. 6.28 % (Bures: 10.32 %) 
were detected by fully decomposable and 0.44 % (Bures: 
1.06 %) by fully PPT witnesses. As expected, fully PPT 
witnesses detect fewer states. 

While the problem can still be tackled numerically for 
six or seven qubits, in recent experiments up to 14 ions 
have been coherently manipulated [Hf. Therefore, we 
study analytical witnesses which can be generalized to 
an arbitrary number of qubits in the following. 

Analytical results — A fully decomposable witness for 
the four-qubit linear cluster state \Cl 4 ) [23| that is ob- 
tained by the SDP of Eq. (J4| is given by 



TF< 



C14 



1 

2 1 



\Gh){Cl A \--(t- 9l ){t-g i ) , (5) 



where g\ = Z1Z2I3I4 and ,94 = li 12^3^4 are two of the 
generators of the cluster state's so-called stabilizer group. 
This witness detects more states than the usual projector 
witness M^proj = |l — \Cl n }(Cl n \, since Wcia is obtained 
from Wp r oj by subtracting a positive operator P + . For 
n qubits, the generators are, after a local basis change, 
gi = XiZ 2 , g% = Zi-iXiZi+i for 1 < i < n and g n = 
Z n -\X n . Then, the n-qubit linear cluster state is defined 
by \Cln)(Cl n \ = 2" n nr=i (! + Thc construction of 
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the four-qubit cluster state witness can be generalized 
to an arbitrary number of qubits (see Appendix E) . For 
seven qubits, e.g., a witness is given by 

w a7 = h- \ci 7 )(ci 7 \ - -L [(i - gi) (i - .94) (i - .97) 

l lb 

+ (1 + gi ) (1 - g 4 ) (1 - g 7 ) 

+ (l-<7i)(l + 54 )(l-£/7) 

+ (t-g 1 ){t-g i )(t+g 7 )] . (6) 

For the case of n qubits, the white noise tolerance is 

Ptoi-(l-2-"+ 1 + (fc + l)2- fe )" 1 ^>l (7) 

where k = L^ipJ- This result is remarkable for sev- 
eral reasons. First, Wci n is the first example of a wit- 
ness for genuine multipartite entanglement so far whose 
white noise tolerance converges to one for an increasing 
number of qubits. Thus, the volume of the largest ball 
inside the biseparable states around the totally mixed 
state approaches zero. A similar scaling behavior of the 
entanglement in the cluster state has been found in Rcf. 
[30j| . Note that, however, they considered full separa- 
bility and not genuine multipartite entanglement. For 
full separability, this scaling behavior is not surprising, 
since it is known that the largest ball of fully separable 
states around the totally mixed states shrinks with in- 
creasing particle number [29| . Moreover, the white noise 
tolerance of Eq. (JT]) corresponds to a required fidelity 
-Freq ~ 1— ptoi ~ k2~ k for large n and therefore decreases 
exponentially fast with growing n. In contrast, the fi- 
delity needed to detect entanglement using Wp ro j equals 
one half, independent of the particle number. Interest- 
ingly, this exponential improvement comes with very low 
experimental costs, since the additional term P + can be 
measured with only one experimental setting. Finally, 
note that our construction induces a similar construction 
for the 2D cluster state. 

Discussion — In this Letter, we presented an easily im- 
plementable criterion for genuine multipartite entangle- 
ment. We demonstrated its high robustness, connected 
it to entanglement measures and provided powerful wit- 
nesses for an arbitrary number of qubits. 

Due to its versatility, the criterion can be used to 
characterize the entanglement in various physical sys- 
tems, e.g. in ground states of spin models undergoing 
a quantum phase transition. Moreover, it is a promising 
tool to study multipartite entanglement in continuous- 
variable systems. Finally, we believe that, as an easy-to- 
use scheme, it will be valuable for the analysis of experi- 
mental data that do not constitute a whole tomography. 

We thank M. Kleinmann, T. Monz, S. Niekamp, A. Os- 
terloh, M. Piani and G. Toth for discussions and acknowl- 
edge support by the FWF (START Prize and SFB FO- 
QUS). 
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APPENDIX 

This appendix consists of five sections. In Appendix 
A, we prove the properties (i) to (iv) of the entanglement 
monotone given in the main text. Appendix B gives an 
example for another superset that can be used to approx- 
imate the set of states which are separable with respect to 
a fixed bipartition. In the third section of the Appendix, 
we prove the optimality of some white noise tolerances 
given in Table I in the main text. Afterwards, we write 
down a fully decomposable witness for the Dicke state 
with two excitations in Appendix D. In the main text, 
this witness was used to illustrate the case in which one 
considers only witnesses that consist of a restricted set of 
observables. Finally, in Appendix E we provide a general 
theory for witnesses of linear cluster states. Examples 
for these witnesses have been given in the main text. 



Appendix A 

In this appendix, we introduce an entanglement mono- 
tone for genuine multipartite entanglement that is moti- 
vated by the notion of fully decomposable entanglement 
witnesses. Let us point out that this defined quantity 
equals the negativity for the bipartite case hence it 
can be considered as an extension of the negativity to the 
multipartite case. 

Proposition 1. For a generic multipartite state g, con- 
sider 



N(g) 



min Tr (nW). 

wew 



W = {W\ VM : 3 < P Ml Q M <1:W = P ai 



(8) 

1- 



where M stands for a possible partition of the subsystems. 
Then N(g) has the following properties: 

• N(g hs ) = for all biseparable states g hs . 

• -/V"[A LO cc(£0] < N(g) for all full LOCC operations. 



N(U loc gU 



loc) 



N(g) for local basis changes U\ c 



• NQ^iPiQi) < ^2iPiN(gi) holds for all convex com- 
binations '^2 i Pigi- 

Proof. The first statement follows directly from the fact 
that any fully decomposable witness can only detect gen- 
uine multipartite entanglement, hence the expectation 
value satisfies Tr(g hs W) > and vanishes for W = 0. 

For the second property we effectively show A[A(g)] < 
N(g) for all trace-preserving, completely positive oper- 
ations A(g) = J2t K *e K l with Y.i K \ K * = 1 that ad- 
mit a fully separable operator-sum representation, which 



means that each operator Ki — Ai <Ei Bi (£) ■ ■ ■ ® Fi has 
a tensor product form. This set of operations defines a 
superset to the set of full LOCC operations, so we verify 
the aboveproperty for an even larger set of possible op- 
erations [2|. Let us point out that for each operation A, 
there exists an adjoint operation A^(Y) = ^ZiK^YKi, 
that satisfies the identity Tr[A(A)F] = Tr[XAt(F)] for 
all linear operators X, Y. The trace-preserving condition 
for A translates to a unital condition for the adjoint map 
At(l) = 1. 

Let us first prove the following statement: For each 
trace-preserving, separable operation A and for any de- 
composable operator W the observable A^ (W) is decom- 
posable as well. Suppose that W = P + Q Tm is an 
appropriate decomposition 0] with respect to a chosen 
partition M. Because of linearity we obtain A^(W^) = 
A^P) + A^(Q Tm ). First, we want to check "normal- 
ization" < A^(P) < 1. Complete positivity provides 
AT(P) > since P > is positive semidefinite itself. If 
one applies the adjoint map to (1 — P) > and employs 
the unital condition one obtains 

A f (l - P) = At(l) - At(P) = l - At(P) > 0, (9) 

hence the upper bound A^ (P) < 1 holds as well. We can 
apply the same argument to A^(Q) if we can fulfil the 
identity A^Q 7 ") = [A i <(Q)] TM with A being completely 
positive and unital. Using the assumed tensor product 
structure of each operator Ki it is straightforward to de- 
duce the Kraus operators of the liner map A satisfying 
this identity. These operators Ki = Ai® Bi® ■ ■ ■ C§) Pj are 
given by A\ = A t if A M and M = A* if A e M, and 
similar relations for all other operators. Let us point out 
that this is the only step where one explicitly needs the 
separable operator structure. 

Via this statement one finally obtains 



N[A(g)} = 



min TrlAfoWl = - min TvlgAUW)} 
wew wew 



< 



min TrlpVFl 
wew 



N[g}. 



(10) 



The inequality holds since A^W) is decomposable again, 
hence the optimization over the complete set of decom- 
posable entanglement witnesses can only produce a lower 
(negative) expectation value. 

Invariance under local basis changes is a direct con- 
sequence of the previous property. Since any local basis 
change U\ oc represents an invertiblc full LOCC operation, 
one obtains 

N(g) > N{U locQ Ul c ) > N[ui c (U loc gUl c )U loc ] - N(g). 
Thus a local basis does not change the value of N(g). 
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The convexity statement 



i 



-mmTr(piW) 
w 



(11) 



follows from linearity of the trace and the fact that if one 
performs independent optimizations then the obtained 
expectation value can only be smaller. This completes 
the proof. □ 



Appendix B 

Here, we would like to give another example of a super- 
set that one can use to approximate the separable states 
of a given bipartition. For the bipartition A\BC, one 
can consider states that possess a symmetric extension 
to k copies of system A, cf. Refs. [J]. This means that 
the given state Qa\bc can be written as the reduced state 
of a multipartite state QAi...A k \BC that is invariant un- 
der all possible permutations of the copied subsystems. 
Every separable state necessarily satisfies this extension 
condition for any number of copies and we denote states 
of this class by Q^^bc' Consequently, we ask whether a 



three-particle state can be decomposed as 



6 - V\Q A \ BC + P^e B \AC + P30 C \AB- 



(12) 



Any biscparable state can be written in this way, hence if 
this expansion fails, genuine multipartite entanglement is 
detected. This approach has appealing properties: With 
increasing number of copies, these supersets converge to 
the set of separable states [J|. Moreover, it is again pos- 
sible to characterize such decompositions using entangle- 
ment witnesses that can be tackled via SDP. These wit- 
nesses are such that for all possible bipartitions M, the 
operator W must be a bipartite entanglement witness for 
the case of k symmetric extensions as given in Ref. 



for a white noise contribution of p « 0.521 the corre- 
sponding state is biscparable. 

Consider the tripartite case. Instead of supersets to 
approximate separable states for bipartition A\BC, we 
now consider a strictly smaller set of states, i.e., any state 
of this subset must necessarily be separable. States of 
this restricted class are denoted as qa\bci an d similar for 
other bipartitions. If a given state has a decomposition 



Q = V\QA\BC + P2QB\AC + P30C\AB, 



(13) 



then it is clearly biscparable. The main difficulty is 
to find operational subset approximations, though re- 
sults have been obtained along this direction recently 
8]. For the thrce-qubit case wc employ an approxima- 
tion that exploits the low dimensionality of the prob- 
lem. Instead of all separable states of a 2 £g> 4 system 
we consider only separable states Qa\bc that are sup- 
ported on the symmetric subspace of system BC . Since 
Sym("HBc) — one effectively has a 2 (X) 3 system for 
which the PPT condition Qa\bjj — becomes necessary 
and sufficient for separability [9j. Similar definitions are 
used for the other bipartitions. The search for an explicit 
decomposition as given by Eq. (fil?]) can be cast into the 
form of a SDP again. 

The solution for g(p) = (1 -p)|W 3 )(W 3 | + pl/8 is as 
follows: The separable states for B\AC and C\AB are the 
same as for the bipartition A\BC but with appropriate 
permutations of the subsystems. All probabilities equal 
Pi = 1/3. Thus, it is left to characterize Qa\bc- The 
corresponding eigenspectrum is given by Ai = A2 = p/8, 
A 3 = 3p/8 and A 4 = 1 — 5p/8 with respect to the eigen- 
states - |000), |x 2 ) = |111), 

| X3 ) = (2a/2|1>®|V+>-|011>)/3, (14) 
|X 4 ) = a|0)®|^ + ) + V / l-a 2 |100), (15) 

with the abbreviation = (|01) + \10})/y/2. The pa- 
rameter a € 1R is determined by the linear condition given 
by Eq. (fl~3|) and will be the root of a non-linear function. 
The only remaining condition is Q T a ^bc — ^ which is sat- 
isfied if p > p with a critical white noise level given by 



Appendix C 

In this appendix, we prove that our criterion is opti- 
mal for three-qubit GHZ and W states and for four-qubit 
GHZ and cluster states (cf. marked states of Table I) . 

First, for the GHZ states mixed with white noise, the 
noise thresholds at which the states become biseparable 
were already derived in Ref. [H| and the values of p to \ for 
our criterion coincide with these values. 

Second, for the three-qubit W state \Ws) we introduce 
a variation of our relaxation idea, in order to show that 



P 



1 

382 



367 - 7lV3 - \/ 2894 - 2988 



0.52102. 



(16) 

At this white noise level one has a w 0.98716. The noise 
level p coincides with the maximal ptoi for our criterion. 

Finally, for the four-qubit cluster state \CU) we used 
the algorithm of Ref. [l0( , that can be used to prove dif- 
ferent forms of multipartite separability. We found that 
a state with a noise contribution of p = 0.616 is bisepa- 
rable, while the witness Wcu of Eq. (6) in the main text 
detects it for p < 8/13 0.6154 as genuine multipartite 
entangled. 
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Appendix D entanglement detection by measuring a restricted set of 

obscrvables is given by 

The fully decomposable witness for the Dicke state 
|Z?2,4) which was used in the main text to demonstrate 



W D = — + a x X m + ai Y m + a 2 Z m + a 3 (X 1 X 2 Y 3 Y 4 + perms) + a 4 (Z 1 Z 2 Y 3 Y 4 + perms) (17) 
16 

+a 4 {Z\Z 2 X 3 Xi + perms)+ a 5 (X1X2I3I4 + perms)+ a 5 (Y1Y2I3I4 + perms)+ a e (ZiZ 2 t 3 ti + perms)] . 

Here, X1X2Y3Y4, + perms is the sum over all distinct permutations of XiX 2 Y 3 Yi. Moreover, ot\ = 0.014, a 2 = —0.095, 
a 3 = 0.0046, a 4 = 0.16, a 5 = -0.14, a 6 = -0.15. 



Appendix E 

Here, we prove the generalized n-qubit version Wcin 01 
the cluster state witnesses given in the main text to be 
fully decomposable witnesses. For four and seven qubits, 
Wcin reduces to forms of Eqs. (6) and (8) of the main 
text, respectively First, let us briefly introduce the no- 
tation that we use. 

Consider an arbitrary graph G = (V, E) that is defined 
by a set V of vertices which correspond to qubits and a 
set E of edges that connect some of these vertices. 
Then, by N(i) we denote the neighborhood of qubit i, i.e., 
the set of all vertices that are connected to i by an edge. 
Moreover, we define N{i) = N(i) U {i}. 
The commuting operators defined by 

g i = X i J] Z k ,i = l,...,n (18) 

fceV(i) 

are the generators of the so-called stabilizer group, i.e., 
{Si, . . . , 52™} = (<7ij • • • ,g n )- In other words, every sta- 
bilizer Si can be written as a product of some (or all) of 
these generators gi. 

The graph state \G) associated to the given graph G is 
then uniquely defined by 

9i \G} = \G), V* = l,...,n. (19) 

One can also associate an orthonormal basis to a graph 
G. The elements |ai . . . a n )o of this so-called graph state 
basis are defined by 

gi\a\ . . .a„) G = (-l) a *|ai . . .a n )o, Vi = l,...,n. (20) 

Consequently, \G) = |0 . . . 0)g?. Also, note that projectors 
on these vectors can be written as 

G \a 1 ...a n ){a 1 ...a n \ G = YY ^ ■ (21) 

i=l 

For better readability, we will drop the subscript G in the 
following and write \ai . . . a n ) q = \a) . Thus, except for a 



single exception that we will explicitly mention, all vec- 
tors in the following have to be understood in the graph 
state basis. However, it is important to note that all par- 
tial transpositions will be done w.r.t the computational 
basis. 

A partition M is a subset of the set of all qubits, while 
we refer to the set made up of partition M and its com- 
plement M as bipartition M\M. 

With the notation clarified, we first prove three lemmas 
to prepare the proof that the operator Wei presented 
in the paper is a fully decomposable witness. The first 
of these lemmas starts from two orthogonal graph state 
basis vectors and shows which kind of partial transpo- 
sitions one can apply to one of these vectors such that 
the two resulting vectors are orthogonal. Lemma [3] sup- 
plies an upper bound on the eigenvalues of any partially 
transposed state in terms of the state's Schmidt coef- 
ficients. Finally, Lemma 2] provides some kind of con- 
struction manual. Given a qubit and starting from a 
graph state basis vector, it tells us how we can construct 
expressions out of this basis vector which are invariant 
under a partial transposition on the given qubit. 
These lemmas will then be used in Proposition[5]to prove 
that the presented operator Wcin is a billy decomposable 
witness. 

Lemma 2. Given a graph G = (V,E) of n qubits and 
an arbitrary bipartition M\M of these qubits. Let \a) 
and \b) be two arbitrary states in the associated graph 
state basis. Let Um be the set of all qubits that lie in the 
same partition as their neighbors, i.e., Um = {i \ •AZ'(i) C 
M V Af(i) C A/}. If there is a qubit i G Um, s.t. b t ^ a*, 
then 

(b\(\a)(a\f'"\b)=0 (22) 

Proof. Let gi,i = l...n be the generators defined by 
Eq. (HI]). Since X T = X, Y T = —Y, Z T = Z and 1 T = 
1, the partial transposition of a product of generators 
only changes the product's sign. Thus, we can describe 
the action of the partial transpose Tm w.r.t partition M 
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on products of generators by Note that the support supp(/) of a Boolean function con- 

tains the bits that the function depends on, i.e., 

/ n \ T " / n \ 

(II,; j ' 1' " (ll'// J • (23) 

where Xi £ {0, 1}. Here, / depends on M and is a 
Boolean function defined by 




supp(/) = {i | 3 f, s.t./(xi, . ..,Xi,...,x n )^ f(xx,- ..,Xi®l,.. .,x n )} . (25) 

Due to the explicit form of the gt, flipping the value of Xi cannot change f(x), Hi £ Um- Therefore, WmHsupp(/) = {}. 
For this reason, we can pull qubits vciIAm out of the partial transposition Tm in the following way (using also Eq. (|2ip). 

{b\(\d){a\f M \b) = Tr 



n 



n 



(-l) Q *. 9i + l 




(26) 



Since 



g.+l g«+l _ |- ne i as ^. expression vanishes if 



2 2 

there is an i £ IAm , such that bi =/= a 



□ 



Lemma 3. Given a state and its Schmidt decompo- 
sition \ip) = YliLi with respect to some bipar- 
tition M\M, where A, > 0, di = dim(Af), d 2 = dim(M) 
and w.l.o.g. d\ < d 2 - Note that, here, \^i)®\vi) is not to 
be understood in a graph state basis. Then, for any state 

\4>), 



<«^><V>|) Tm |^> < maxA| 



(27) 



Proof. Writing down (ip |) M in the basis {\/J>i) 65 
\ v i)}i=i...d 1 ,j=i...d 2 ) one obtains a matrix with two dif- 
ferent kinds of submatrices. First, a diagonal one with 
diagonal elements Af or zero. Second, anti-diagonal sub- 
matrices of the form 







XiXj 







(28) 



Thus, the eigenvalues of the total matrix are 
{±AjAj, Af , 0} and the maximum of these eigenvalues has 
the form Af. 

□ 

For the following lemma, we need to keep in mind that 
the application of the Pauli operator Z k to a graph state 



basis vector creates a bit flip on bit k, i.e., 

Z k \a) = |ai . . . a k -i flfc © 1 a k +i ■ . . a n ) . (29) 

Lemma 4. Given a graph G. Then, in the associated 
graph state basis, 

\d)(d\ + \b)(b\) Tk ^\d)(a\ + \b)(b\, (30) 



if \b) is obtained from \a) in one of the two following 
ways. 

(i) \b) = Z k \S). 

(ii) \b) = n ieA f(fc) Z i\ E )- 
Proof. (i) With Eq. ((2TJ) and Eq. flU]), we have 

\d){d\ + \b){b\ 
'g k + t , -g k + l\ A (-l) a 'g, + t 



2 2 J 2 

i=i 

i^k 



n 



n 

i=l 
i^k 



(31) 



Since g k cancels in Eq. (|3"Tj) , the explicit form of the 
generators gi implies that, in this equation, there is 
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no Y on qubit k. Since Y is the only Pauli matrix 
that changes under partial transposition, |a)(a| + 



Zk\a){a\Zk is invariant under Tfc. 



(ii) Using again Eqs. (HQ) and ^ yields 



m+M\= n 



(-l) a * 9i + l { (-l) a ' +1 ff t + l j (-l) a ", 9fc + l (-l) a '.g t + 1 

»GA/"(k) / i£AA(fc) 



(32) 



The expression in the first brackets can be simplified to 
j-j (-1)^ + 1 , tt (-l^ft + l 



ieN(k) 



n 

ieJV(fc) 



2 ~ n \ e n 9?+ e (-i) a " +E * x - n 9i 

\x£{0,l} N i&N{k) xe{0.1} N i&N{k) 

^ m-f n sf'l > ( 33 ) 

Veven iejV(fc) 



where we defined ./V = |A/"(fc)| and Xieven * s * ne sum over ai l & G {0, 1} W for which 5Z i=1 £j is 
Combining Eqs. fl32j| and fl35|) yields 

\a ) (a\+\b)w=2- N (j2(- i r n sr)(-ir^ n ( ~ 1)n ' 2 g<+1 

\ cven ieW(fc) / i£AA(fc) 

\even ieAA(fe) / i£A?(fc) 



even. 



(34) 



r 



In Eq. (|34p , the right-hand side is invariant under 
Tfc. To see this, consider the first part of the sum, 
in which all terms contain and an even number 
of generators in Af(k). Together with the form of 
the generators [see Eq. (fTB)) ] , this implies that there 
is no Y on qubit k. 

The second part of the right-hand side of Eq. ([34]) 
does not have a Y on qubit k either, since it does 
not contain g^. Consequently, Eq. AMI) is invariant 
under Tfc. 

□ 

Let us now return to the main proof in which we will 
need Lemmas [21 [3] and [4] 

Proposition 5. Let \Cl n ){Cl n \ be an n-qubit linear clus- 
ter state with n > 3 and gi the generators of its stabilizer 
group. Let B — {Pi} be a subset of all qubits such that 
fif(Pi) r\Af(j)j) = {} for i + j. We define m = \B\. Let 
^ s be the sum over all vectors s of length m with ele- 
ments Si = ±1 that contain at least two elements that 
are —1, i.e., Xa=i Si < m — A. Then, 



1 



W cln = -l-\Cl n ){Cl n \ 



En 



Sigi 



(35) 



is a fully decomposable witness for \Cl n )(Cl n \, i.e., for 
any strict subset M of all qubits, it can be written in the 
form Wet = Pm + Qlf , where P M > 0, Q M > 0. 



For 



the sake 

Sjgi + 1 



of brevity, we define P+ = 
Note that P + is a sum of all projec- 



ts ILee 2 

tors onto graph state basis vectors that contain at least 
two excitations in B, i.e., two bits that equal one. For 
example, if we choose B — {1,4, 7, . . .} for the sake of 
illustration, 

P+ = (IOX1X2IX3X4I . . . ) (0x1X21x3X4! . . . j 



IIX1X2OX3X4I . 
IIX1X2IX3X4O . 
IIX1X2IX3X4I . 



. )(lXiX20x3X4l . . . 
. )(lXiX2lX3X40 . . . 
. )(lXiX 2 lx 3 X4l . . . 



+ ...) 



(36) 



This also illustrates why the presented construction does 
not work for linear cluster states that consist of three or 
less qubits. In such a case, B could only have one element, 
since the elements' neighborhoods cannot overlap. Then, 
however, one cannot fulfill the condition that there must 
always be two bits in B that equal one. Note that sets 
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B = {1,4} and B = {1,4,7} for four and seven qubits, 
respectively, result in the witnesses Wen and Wen 01 
Eqs. (6) and (8) of the main text. Moreover, these choices 
of B lead to the noise tolerance of Eq. (9) of the main 
text. 



In other words, r equals the number of bits Pi € B which 
do not lie at the border of M (or M), i.e. they do not 
obey 



Af(A-) C M or AT(Pi) C M. 



(39) 



Proof. First, we will provide a way to construct an ap- 
propriate operator Pm > for every M. The main trick 
is to do this in such a way that this operator is invariant 
under certain partial transpositions. Second, we will use 



{W-Pi 



M) 



> for ev- 



these Pm to show that Qm 
ery M. 

Note that we order the qubits Pi in a canonical way such 
that Pi < pi+\. To simplify notation, we define bitstring 
wi for every Af(Pi) (of \Af(Pi)\ bits length) in such a way 
that its bits define to which partition the corrcspoding 
qubits belongs. w*i = (oj-i, uio, w+i), w j G {0,1}, means 
that qubit Pi + j does not belong to M if uij = and it 



For Pi being the first qubit, we 
for Pi being the last qubit, we 



belongs to M if u)j = 1 . 
consider Wi = (cjo,^+i) 
use Wi = (cj-i,cj ). 

Given M , we proceed in the following way to construct 
Pm- 



1. Start with P, 



(0) 



M 



\Cl n )(Cl n \ = 0...0 0...0. 



2. Loop through the values ofz = l,2,...,m and up- 

( i) 

date P M until you reach i > m. Then proceed 

with step 3. The update of is carried out in 
the following way. 



(i) If = 2: 

I If m = 01 or Wi = 10 

P M ~ 1] + Z k P M ~ 1] Z k , where k e 7V(ft) 
II If (Z;; = 11 or m 

(ii) If \ti(Pi)\ = 3: 



then P$ = 



00, let P« = P^~ 



I If ^ = 110 or uJj = 001, then P 



(i-i) 



II If m, 



(i) 
M 



M 



J /3i + l- 



= 010 or m = 101, then P 



(i) 



(i-i) 



III If Wi 



100 or w. 



lZ/3, + lP M 1) Zfj i 



011. then P 



P 



(i-i) 



M 
(i) 



P 



IV In all other cases, let P^ — -i M 

Let r be the number of times that one of the cases 
2.(i).I, 2.(ii).I, 2.(ii).II or 2.(ii).III occurred, i.e., the 

number of steps in which P M changed. 
If r < 1, define 



M 







(37) 



Let t be the value of i for which PjV was changed 
the last time, i.e., P$ 
define 



P M ] V i > t. If r > 1, 



M 



P. 



(t-1) 



M 



\Cl n )(Cl r 



(38) 



Thus, for r = 0, all Pi £ B obey this property, for r = 1, 
there is one exception. 

Note that the operator P M constructed via the given al- 
gorithm is either zero or a sum of one-dimensional pro- 
jectors onto basis states, i.e., 



P 



M 



(40) 



This can be seen by the fact that P^f — \Cl n )(Cl n \ = 
|0 . . . 0) (0 . . . 0|, the application of Z only flips a bit and 
finally \Cl n )(Cl n \ is subtracted again. Let us illustrate 
the algorithm by a concrete example. 
Example of the algorithm: Consider an eight-qubit 
cluster state and M = {2, 3, 5}. We choose B = {1, 4, 7}. 
Then, the algorithm proceeds as follows. 



100000000) (000000001 



P$ - \Cl 8 )(Cl 8 \ 

Step i = 1, i.e., Pi = 1: Since w\ = 01, case 2. (i) .1 
applies and 



P (i) _ 

M — 



100000000) (000000001 

+ Z 2 |00000000)(00000000|Z 2 

100000000) (000000001 

+ |01000000) (010000001 . 



(41) 



Step i = 2, i.e., 
2. (ii) .II applies and 



4: Since w 2 = 101, case 



P 



(2) _ p(l) , 7 7 p(l) 7 7 

= 100000000) (000000001 
+ |01000000) (010000001 
+ 100101000) (00101000) 
+ |01101000)(01101000| 



(42) 



Step i = 3, i.e., p 3 = 7: We have w 3 = 000. Thus, 
case 2.(ii).IV applies and Pff = Pff- 

Asi = 4>m = 3, we abort the loop. P^ changed 

(i) 

in two steps, i.e., r = 2. Moreover, since 

P M did 

not change in the last step, i.e., the third step, t = 
2. Then, 



M 



P. 



(t-1) 



M 



Ch)(Cl 7 \ 

|00000000)(00000000| 
= |01000000)(01000000| 



= p, 



M 



(43) 
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Let us now return to the general case and understand the 
properties of the operator Pm for an arbitrary M. The 
construction uses Lemma 0] to ensure that, in every step, 
either 



or 



P 



M 



Tm- 



P 



M 



(44a) 



(44b) 



hold, where we defined M k = M nAf(j3 k ). 
To see that Eqs. (|44|) hold, it is best to consider an ex- 
ample: Assume that io, = 010. Then, case 2.(ii).II of the 
algorithm applies. As seen above, P M ~^ nas the form 



P 



(i-i) 

M 



Eb\b)(b\. Thus, 



M — 



+ Zi-iZ i+ i\b){b\Z i+ \Zi. 



(45) 



Then, Lemma 2.(ii) implies that Eq. ()44b[) holds. 
Eqs. ([44|) hold in every step, i.e., for i = j and for 
i = k, where j ^ k. According to the premise of non- 
overlapping neighborhoods of the qubits in B, we have 
Af(/3j) n M{(3k) = {}■ Therefore, the partial transposi- 
tions in Eqs. (|44|) for i = j always affect qubits different 
from the ones that are affected by the partial transpo- 
sitions for i = k. For this reason, Eqs. (|44p for P M 
hold with respect to to every value of k, except for k = t. 
More precisely, 



P 



(t-i) 

M 



= P 



At 



-l)\ T JV(/3 fc ) 



M 



P 



(t-1) 



M 



(46a) 



(46b) 



is true for every k ^ t. We will use this important prop- 
erty later. 

Let us proceed with the proof. Since Pm is zero or has 
the form of Eq. ({40]). we know that P M > 0. Thus, it 
remains to show that (W C i - Pm) Tm = Qm > 0, V M. 
Operators that are diagonal in a graph state basis can 
be written in the form dV i where the sum runs 

over binary strings x that depend on the operator. Since 
any partial transposition can at most introduce minus 
signs in some terms of this sum, such operators remain 
diagonal under any partial transposition. Therefore, it is 
enough to prove that 



(Wei -P M f M \k) >0, VM,\k) 



(47) 



Note that P+ = J2sY[ieB Sl9 '^^ is invariant under any 
partial transposition, since it does not contain any Y op- 
erators due to the form of the generators and the fact that 
no two qubits in B have a common neighbor. Therefore, 



together with the explicit form of the witness given in 
Eq. (f35|) . we can rewrite Eq. (|47|) as 



l-l(k\P + \k)-(k\ (\Cl n )(Cl n \+P M f M \k) > 0, VA/,|fc). 

(48) 

In order to prove this, we distinguish two different cases: 
1. 



\P + \k) ^ ^ (k\P+\k) = I 



Note that the equivalence can best be see in 
Eq. (|36|) . Also, Eq. (|36)l implies that there must 
be at least two bits of k, say ki a ,kj a € B, with 
io 7^ jo 7 such that ki — kj Q = 1. 
In the case P M = 0, Eq. ((48]) and (k\P + \k) = 1 are 
equivalent to 



(k\(\Cl n )(Cl n \f M \k) >0VM,|fc) 



(49) 



To see that the left-hand side always vanishes, one 
uses that Pm = is equivalent to r < 1, i.e., cases 
2.(i).I or 2.(ii).I or 2.(ii).II or 2.(ii).III in the al- 
gorithm occurred at most once. This means that 
N{Pi) C M or Af(Pi) C M holds for all qubits 
G B with at most one exception, namely fit- 
With ki = kj a = 1, Lemma 1 can be applied to see 
that Eq. (|49|) vanishes. 

In the case Pm ^ 0, Eq. (|4"5)l can be simplified using 
(k\P+\k) = 1 to 



(\Cl n ){Cl n \+P M y M \k) >0VM,|fc) 



,(*-!) 
M 



\k) > 0VM, |fc) 



(50) 



Here, the definition of Pm, Eq. (|38|) . has been used 



Now, Pm and therefore P- 



(t-i) 



consists of a sum 



of projectors onto graph basis states 



(see 



Eq. (|3D]l). Since the algorithm starts with Pff = 
|0 ... 0) (0 ... 1 and never flips any bits on the 
qubits Pi G B, these states \a) obey ap i =0, Wi = 
l,...,m. Also, depending on whether io = t or 
jo = t, Eqs. (|46| can be applied to whichever of 
these two qubits is different from t. The invariance 
given by Eqs. (|46|) then implies that this qubit can 
be treated as if it lay on the border of M (or AI) 
and therefore Lemma [2] yields 



(*l [P. 



)(*-!) 
M 



\k) =0 



(51) 



and therefore Eq. ([50|) holds. 



|P+|fc) = 



To show that Eq. (|48|) holds, we need to prove that 
(k\(\Cl n )(Cl n \+P M f M \k)<\- (52) 
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In the case Pm ^= 0, Pm is given by Eq. (|40j) and 
Eq. (|52|) is equivalent to 

(\ T M 
\Cl n )(Cl n \+^2\a)(d\\ \k)<\. (53) 

Note that \Cl n )(Cl n \+J2a = p m~ 1] consists 
of 2 r terms, as one starts with one term and dou- 
bles this number (r — 1) times to obtain P M 1 \ 
Therefore, it is enough to show that 

(k\(\Cl n )(Cl n \f M \k) <2- r (54) 

and 

(k\(\a)(a\f M \k) < 2" r V|a) (55) 

holds. We will do this by using Lemma [3] How- 
ever, since the vectors \a) are basis vectors of the 
graph state basis of \Cl n ), \a) and \Cl n ) are LU- 
equivalent. Therefore, they have the same Schmidt 
coefficients and Lemma [3] results in the same upper 
bounds. For this reason, it suffices to show only 
one of these upper bounds, namely 

(k\(\Cl n )(Cl n \f M \k) <2~ r . (56) 

The idea behind the subsequent proof of Eq. (f56|) 

(i) 

is that, every time P M changes in our algorithm, 
r is increased, but at the same time, a Bell pair 
which has Schmidt coefficients {-^= , ^= } is created 
which results in a smaller upper bound according 
to Lemma [3] 

To easily calculate the Schmidt coefficients of \Cl n ) 
with respect to bipartition M\M, we transform 
\Cl n ) via local transformations into disconnected 
Bell pairs whose qubits arc in different partitions. 
In order not to change the Schmidt coefficients dur- 
ing this transformation, we apply the following two 
operations to the graph. 

• Local Complementation (LC): Pick a qubit j 
and invert its neighborhood graph. In other 
words, if between two neighbors of j there is 
an edge, delete it. If between two neighbors 
of j there is no edge, add an edge between 
them. Since this operation corresponds to the 
application of local unitaries, the Schmidt co- 
efficients remain unchanged. 

• Controlled-Z (CZ): Apply a controlled-Z op- 
eration to two qubits in the same partition. If 
they are connected by an edge, this will delete 
the edge. If they are not connected, this will 
create an edge. As this unitary operation is 
only applied on qubits in the same partition, 
it does not change the Schmidt coefficients. 




FIG. 2: For every case in which Pj} changes in the algorithm, 
we explicitly demonstrate that it is possible to create a Bell 
pair via a sequence of local transformations. Note that the 
cases Wk = 100 and Wk = Oil follow from Wk = 110 or Wk = 
001 due to symmetry reasons. 

Since these transformations can be applied to each 
set Af(f3j ) individually to create a Bell pair between 
the two partitions (see Fig. [2J it is possible to end 
up with distinct Bell pairs between the two parti- 
tions. 

Moreover, the number q of Bell pairs that one can 
create by transforming the qubits in each N{fij) 

separately obeys q > r, as every time that P}J 
is changed in the algorithm (including step i = t) 
corresponds to a configuration shown in Fig. [2] 
Furthermore, each of these configurations allow for 
the creation of a disconnected Bell pair. 
Since the Schmidt coefficients of a single Bell pair 
are , ^= }, Lemma [3] yields 

(k\(\Cl n )(Cl n \f M \k) < 2-9, (57) 
With q > r, this implies that Eq. ([5rJ|) holds. 
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In the case Pm = 0, we need to show that 

{k\{\Cl n ){Cl n \f M \k)<\. (58) 

Note that Pm — is equivalent to r < 1. On 
one hand, r = means that for each i = 1, . . . , n, 
either Af(i) C M or 7V(i) C M holds. Since 
A/ 7^ {1, . . . ,n}, deleting all edges within M and 
M will lead to at least one Bell pair, which has 
Schmidt coefficients , -j= }. With Lemma El 
Eq. (JSHl) holds. 

On the other hand, for r = 1, Eq. (|58|) also holds, 
since due to Fig. [5] the Schmidt coefficients arc 
given by at least one Bell pair. Therefore, the max- 
imal coefficient is i. 
This finishes the proof of Proposition [SJ 

□ 
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